Pharmaceuticals 2014, 7, 695-722; doi:10.3390/ph7060695 



OPEN ACCESS 



pharmaceuticals 

ISSN 1424-8247 

www.mdpi.com/journal/pharmaceuticals 

Article 

Conformational Analysis, Molecular Structure and Solid State 
Simulation of the Antiviral Drug Acyclovir (Zovirax) Using 
Density Functional Theory Methods 

Margarita Clara Alvarez-Ros and Mauricio Alcolea Palafox * 

Departamento de Quimica-Fisica I. Facultad de Ciencias Quimicas, Universidad Complutense, 
Ciudad Universitaria, Madrid 28040, Spain; E-Mail: maralva32@yahoo.es 

* Author to whom correspondence should be addressed; E-Mail: alcolea@ucm.es; 
Tel.: +34-913-944-272. 

Received: 16 January 2014; in revised form: 22 May 2014 / Accepted: 26 May 2014 / 
Published: 6 June 2014 



Abstract: The five tautomers of the drug acyclovir (ACV) were determined and optimised 
at the MP2 and B3LYP quantum chemical levels of theory. The stability of the tautomers 
was correlated with different parameters. On the most stable tautomer Nl was carried out a 
comprehensive conformational analysis, and the whole conformational parameters (R, [3, 
O, cpi, 92, cp3, 94, 95) were studied as well as the NBO Natural atomic charges. The 
calculations were carried out with full relaxation of all geometrical parameters. The search 
located at least 78 stable structures within 8.5 kcal/mol electronic energy range of the 
global minimum, and classified in two groups according to the positive or negative value 
of the torsional angle 91. In the nitrogen atoms and in the 02' and 05' oxygen atoms of the 
most stable conformer appear a higher reactivity than in the natural nucleoside 
deoxyguanosine. The solid state was simulated through a dimer and tetramer forms and the 
structural parameters were compared with the X-ray crystal data available. Several general 
conclusions were emphasized. 

Keywords: acyclovir; 9-(2-hydroxyethoxymethyl) guanine; tautomer; antiviral drug; 
conformational analysis; Zovirax 
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1. Introduction 

Acyclovir (ACV, 9-(2-hydroxyethoxymethyl) guanine, Zovirax [1,2], Figure 1), is a guanine 
derivative possessing antiviral activity and commonly used in the treatment of herpes. It is a potent 
antiviral agent that is used as a highly specific inhibitor of herpes viruses (HSV) types 1 and 2 [3-7]. 
A series of new guanine base modified tricyclic analogues of ACV and ganciclovir were evaluated for 
activity against herpes simplex virus type 1 and 2, showing similar antiherpetic potency as the parent 
compounds ACV and ganciclovir [8]. The antiherpetic activity was found to be strongly dependent on 
the nature and esteric demands of the substituents in the 6 and/or 7 positions [9]. 

Figure 1. Molecular structure and definition of the torsional angles in tautomer Nl of acyclovir. 



Most persons who are infected with human immunodeficiency virus type (HIV-1) are also infected 
with herpes simplex type2 (HSV-2), which is associated with increased plasma and genital levels of 
HrV-1. Thus, it has been reported that ACV inhibits HIV upon human herpes virus (HHV) coinfection 
in tissue cultures. This activity was found to be correlated with the phosphorylation of the parent drug 
to the monophosphate form mediated by HHV-encoded kinase [10]. Recent studies show that ACV 
decreases the HIV-1 RN and suppress both viruses in coinfected tissues [11], with great impact on 
HIV and HSV-2 seropositive patients [12,13]. The treatment of ACV- resistant HSV infections [7,14] 
and tests for release it at oral therapeutic levels have been also reported [15]. 

As it is shown above, ACV has been extensively studied from the pharmaceutical and medical point 
of view. However, few studies appear on its molecular structure, only at low level [16], and to our 
knowledge, there is no data on its conformational characteristics by Density Functional Methods 
(DFT) or ab initio quantum chemical methods. Thus, it is one of the tasks of the present manuscript. 
Conformations of some derivatives of ACV with biological activities have been reported [17], as well 
as studies of complexes of ACV with several metals [18,19]. 

The tautomeric study of DNA structural components has a great interest today with numerous 
research publications [20,21]. The tautomeric equilibrium in ACV between the keto and enol forms 
has been observed from the UV/VIS spectra [22]. This equilibrium depends on the polarity of the 
solvent, and therefore, in water solution the keto form prevails, while in methylene chlorides it is the 
enol one. Thus, another goal of the present manuscript is to study the possible tautomers of ACV, and 
the determination of their % populations at room temperature and at 273.15 K. The stability of these 
tautomers and its dependence of different parameters is another point analyzed here. 
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Finally, from our understanding, would be interesting to calculate the different conformational 
possibilities of ACV and compare the results with the natural nucleoside deoxyguanosine (dG). This is 
the last goal of the present manuscript. An accurate knowledge of the different conformers of ACV, its 
charge distribution, inter- and intra-molecular interactions, solid structure, and flexibility would be an 
important help for the interpretation of drug-target interactions, as well as to design new antivirus. 

For this reason, the conformers of natural and analogues nucleosides have been analyzed by 
different authors. Now, an extensive theoretical study of the conformational preferences in ACV has 
been carried out with full relaxation of all geometric parameters, in an attempt to gain insights into 
molecular features responsible for activity. We will attempt to determine herein, if the various 
geometric features in ACV are correlated or interact with one another. We are also interested in 
whether alternative forms of hydrogen bonding make significant contributions to the conformational 
behavior of ACV. 

2. Computational Details 

Calculations were carried out by using the Becke exchange functional (B) [23], Becke's 
three-parameter exchange functional (B3) [24], Handy's OPTX modification of Becke's exchange 
functional (O) [25,26] and the extended (X3) [27], in combination with the correlation functionals of 
Lee, Yang, and Parr (LYP) [28], and Perdew and Wang's 1991 (PW91) [29]. B3LYP Density 
Functional method (DFT) is the most used today, and for this reason the majority of the calculations 
were carried out with it. 

All the methods appear implemented in the GAUSSIAN 03 program package [30]. The UNIX 
version with standard parameters of this package was used in the alpha computer of the Computational 
Centre from University Complutense of Madrid, in which all quantum chemical computations were 
performed. DFT methods provide adequate compromise between the desired accuracy and the heavy 
demands put on computer time and power. Different studies have shown that the data obtained with 
DFT methods are in good agreement with those obtained by expensive computational methods as 
MP2 [31-35]. Also they give good results in many studies on nucleosides and on drug design [36,37] 
and they predict vibrational frequencies of DNA bases better than HF and MP2 [38^12]. Moreover, 
because of the high size of ACV, MP2 calculations were only possible with the 6-31G(d) basis set for 
memory computer problems. However, we used the 6-3 lG(d,p) basis to calculate single point energies. 

The 3D Potential Energy Surface (PES) of this molecule was determined by rotation of the 
exocyclic and endocyclic torsional angles (pi (equivalent to % angle in nucleosides [43]), q>2, q>3, cp4 
and cp5. These dihedral angles were simultaneously hold fixed at values varying between 0° and 360° in 
steps of 60° in a first study. All other geometrical parameters were relaxed during these optimizations. 
78 optimized geometries were obtained in this step by minimizing the energy with respect to all 
geometrical parameters without imposing molecular symmetry constraints. Berny optimizations 
under the TIGHT criterion were used. Atomic charges were determined with the Natural NBO 
procedure [44]. 

Harmonic wavenumber computations were carried out at the same level of the respective 
optimization process and by the analytic evaluation of the second derivative of the energy with respect 
to nuclear displacement. Vibrational frequencies calculations were performed on all optimized 
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conformers to confirm that they corresponded to local minima. All optimized structures showed only 
positive harmonic vibrations (local energy minima). Relative energies were obtained by adding 
zero-point vibrational energies (ZPEs) to the total energy. For the calculation of the ZPEs, the 
frequencies were retained unsealed. The AG values were sums of electronic and thermal Free Energies. 
The conformational equilibrium at 298.15 K was evaluated by means of the Boltzmann distribution 
formula exp(-AG/kT), where AG is the relative Gibbs energy. 

3. Results and Discussion 

ACV has five possible tautomers (Figure 2), that were fully optimized at different levels of 
computation, Table 1. The most stable one corresponds to Nl, and thus we have focused the study only 
in this tautomer. The remaining forms were left for future research. In the last two column of this 
Table is shown the % population of the different tautomers at 298.15 K and 273.15 K. At room 
temperature the largest population corresponds to tautomer Nl (48.1%). The second population is due 
to OHC tautomer (37.7%) and the third one is to OHT (14.3%). Tautomers N3 and N7 have very little 
population, less than 0.05%. 



Figure 2. Five tautomers described in acyclovir. 




OHC OHT 
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Table 1. The optimum stable tautomers calculated in acyclovir molecule at the levels: B3LYP/6-31 1++G(3df,pd) (values in bold); 
MP2/6-31G(d,p) (values in key); B3LYP/6-31G(d,p) (values in normal); 03LYP/6-31G(d,p) (values in brackets); B3PW9 1/6-3 lG(d,p) 
(values in parenthesis); B3LYP/cc-pVDZ (values in quotation marks) and B972/6-31G(d,p) (values in italic type). Torsional angles in degrees, 
dipole moments in debyes, distance R in A, energy increments in kcal-moi -1 and population (%) at 298.15 K (P298.15) and at 273.15 K (P273.15). 



Tautomers 


R 


P 


O 


<Pi 


q>2 


93 


q>4 


<Ps 


H 


AE 


AG 


1*298.15 


P273.15 


Nl 


3.963 


72.5 


54.3 


-75.3 


146.4 


-88.5 


69.8 


-69.3 


5.739 


0 a 










{3.892} 


{73.5} 


{51.3} 


{-71.9} 


{150.3} 


{-94.6} 


{69.1} 


{-75.4} 


{6.207} 


{0} b 










3.925 


72.3 


53.7 


-74.1 


144.7 


-93.4 


71.1 


-73.3 


5.375 


r\ C 

0 


0 


48.1 


49.2 




[3.956] 


[72.7] 


[53.9] 


[-75.4] 


[144.5] 


[-94.0] 


[71.5] 


[-72.1] 


[5.369] 


[or 


[0.638] f 








(3.891) 


(72.7) 


(53.2) 


(-75.1) 


(145.2) 


(-92.3) 


(70.3) 


(-74.1) 


(5.492) 


(0) g 


(0) 








"3.922" 


"72.0" 


"54.0" 


"-73.9" 


"143.2" 


"-92. 0" 


"72.3" 


"-74.4" 


"5.34" 


"0" ' 


"0" j 








3.913 


72.8 


53.3 


-75.0 


145.3 


-93.4 


70.8 


-73.6 


5.502 


o k 


o l 






N3 


{3.672} 


{77.0} 


{53.5} 


{-45.9} 


{149.1} 


{-92.5} 


{60.2} 


{-92.6} 


{12.220} 


{25.457} 










3.877 


80.5 


49.5 


-48.4 


155.0 


-93.9 


61.1 


-88. 


11.183 


25.529 


24.811 


<0.05 


<0.05 


N7 


{3.741} 


{73.0} 


{52.0} 


{-76.4} 


{140.5} 


{-80.8} 


{72.4} 


{-92.7} 


{9.841} 


{19.302} 










3.745 


72.8 


53.3 


-79.9 


136.9 


-78.5 


71.8 


-91.5 


7.785 


17.516 


17.689 


<0.05 


<0.05 


OHC 


{3.826} 


{73.5} 


{51.7} 


{-72.3} 


{146.5} 


{-89.2} 


{70.4} 


{-82.6} 


{2.471} 


{0.603} 










3.839 


72.9 


53.4 


-75.9 


141.9 


-87.0 


71.1 


-81.5 


2.100 


0.145 


0.284 


37.7 


37.7 


OHT 


{3.792} 


{73.5} 


{52.0} 


{-72.4} 


{145.0} 


{-87.2} 


{70.6} 


{-86.1} 


{3.437} 


{1.131} 










3.806 


72.9 


53.6 


-76.2 


140.7 


-85.0 


71.0 


-84.9 


2.105 


0.719 


0.830 


14.3 


13.1 



a AE = 0 = -811.214731 a.u. at B3LYP/6-311++G(3df,pd) level; b AE = 0 = -808.643998 a.u. at MP2/6-31G(d,p) level; c AE= 0 = -810.723396 a.u. at B3LYP/6-31G(d,p) level: 
d AG = 0 = -810.764701 a.u.; c AE = 0 = -810.449499 a.u. at 03LYP/6-31G(d,p) level; f AG = 0 = -810.490782 a.u.; g AE = 0 = -810.419555 a.u. at B3PW91/6-31G(d,p) level: 
h AG = 0 = -810.460766 a.u.; ' AE = 0 = -810.760402 a.u. at B3LYP/cc-pVDZ level; j AG = 0 = -810.801662 a.u.; k AE = 0 = -810.444324 a.u. at B972/6-31G(d,p) level; 
'AG = 0 = -810.485454 a.u. 
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An analysis of the relative energies of these tautomers shows that they can be related to the dipole 
moment (u) and to the torsional angle cpi, Figure 3. Thus, the least stable tautomer in the isolated state 
(N7) has the highest \i, i.e., it is the most stable in water solution. Other relations can be observed 
between AE and cpi. It is noted that tautomer N7 has the highest negative value of cpi angle. 

Figure 3. Relationship/tendency observed between the relative electronic energy 
AE + ZPE correction of the different tautomers vs. (a) the dipole moment u. and (b) the 
exocyclic torsional angle cpi. 
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3.1. Definition of the Conformational Angles 

The atomic description of tautomer Nl of ACV, as well as the most important exocyclic and 
endocyclic torsional angles, is defined in Figure 1. The conformation can be characterized by the 
following five torsional angles: (i) the torsional angle, (pi(C4-N9-Cl'-02'), which determines the two 
orientations of the base relative to the opened chain, denoted in the present manuscript as conformers 
A and B; (ii) The exocyclic torsional angle q> 2 (N9-Cl'-02*-C3'); (iii) q> 3 (Cl'-02'-C3'-C4*) that 
determine the different folds of the chain; (iv) q>4 (02'-C3'-C4'-05') describing the 05'H orientation; 
and (v) the torsional angle cps (C3'-C4'-05'-H5') that define the orientation of the hydroxyl hydrogen 
H5'. This cps angle has some similarities with the [3 angle of the nucleosides [43]. 

Considering the structure of the chain, that likes the structure of the sugar in the nucleosides, and in 
accordance to previous works [45], other three structural parameters were defined to fix the chain 
position respect to the plane of the nucleobase: (vi) The vector R (N9---05') which determines the 
distance of the OH group relative to the base; (vii) The angle [3 (C4-N9---05') which defines the angle 
of the OH group relative to the base plane; and (viii) the angle O (CT-N9- • -05') which also determines 
the position of the OH group. 

3.2. Conformers and Energetics 

An extensive conformational study of tautomer Nl was carried out through a rotation of the 
exocyclic (pi, (p 2 , (p3, q>4 and q> 5 torsional angles. A detailed collection of the most important 
conformational parameters of these optimized forms is included in Table 2. The conformers were 
classified according to the two ranges of rotation of q>i: conformers A with the (pi values negative and 
conformers B with (pi positive. 

Two energy criteria were considered for each conformer: the electronic energy AE + ZPE 
correction, and the Gibbs energy AG. For the numbering of the conformers in each range of rotation of 
(pi was followed the AE + ZPE criterion. Calculations at different levels, as well as single point 
calculations at the MP2/6-31G(d,p)//B3LYP/6-31G(d,p) level were carried out to confirm the stability 
of the main conformers, Table 2 and Table 3. In general, the stability order remains, although several 
changes are observed. Thus, conformer B2 appears as the most stable one instead of Al predicted by 
the DFT methods. The % population of the different conformers at 298.15 K (P298.15) and at 273.15 K 
(P273.15) was calculated. It indicates that only confomers Al (41.5%), Bl (36.1%), B2 (9.6%) and A2 
(8.7%) have importance. The population of the remaining conformers is lower than 1% and thus they 
are not of interest. The temperature effect is not significant. 

By the same methodology, the relative energies of the different conformers have been determined 
in related nucleosides [34,46-48]. The global minimum calculated in these molecules by 
MP2/6-31G(d,p) was in accordance to that found by B3LYP/6-31G(d,p). Thus, our results by B3LYP 
can be considered acceptable. 
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Table 2. The 78 optimum stable conformers calculated in tautomer Nl of acyclovir molecule at the levels: B3LYP/6-31 1++G(3df,pd) (values 
in bold); MP2/6-31G(d) (values in keys); B3LYP/6-31G(d,p) (values in normal type); 03LYP/6-31G(d,p) (values in brackets); B3PW91/6- 
31G(d,p) (values in parenthesis); B3LYP/cc-pVDZ (values in quotation marks) and B972/6-31G(d,p) (values in italic type) level. Distance R 
in A, torsional angles in degrees, dipole moments in debyes, and energy increments in kcal-mol -1 . 



Conformers 


R 


P a 




<Pi 


<P2 


<Ps 


<N 


<Ps 


t 1 


AE 


AG 




3.963 


72.5 


54.3 


-75.3 


146.4 


-88.5 


69.8 


-69.3 


5.739 


0 






{2.859} 


{73.8} 


{50.9} 


{-71.8} 


{151.8} 


{-96.3} 


{68.8} 


{-74.2} 


{6.182} 


{0.064} 






3.925 


72.3 


53.7 


-74.1 


144.7 


-93.4 


71.1 


-73.3 


5.375 


0 c 


0 d 


Al 


[3.956] 


[72.7] 


[53.9] 


[-75.4] 


[144.5] 


[-94.0] 


[71.5] 


[-72.1] 


[5.369] 


tor 


[0.638] 




(3.891) 


(72.7) 


(53.2) 


(-75.1) 


(145.2) 


(-92.3) 


(70.3) 


(-74.1) 


(5.492) 


(0) f 


(0) g 




3.922 


"72.0" 


"54.0" 


"-73.9" 


"143.2" 


"-92. 0" 


"72.3" 


"-74.4" 


"5.343" 


"0" h 


"0" 1 




3.913 


72.8 


53.3 


-75.0 


145.3 


-93.4 


70.8 


-73.6 


5.502 


0' 


0 




3.774 


76.8 


52.5 


-97.7 


65.9 


74.6 


-67.0 


-34.1 


5.367 


1.258 






{3.727} 


{77.7} 


{50.7} 


{-96.9} 


{66.0} 


{76.2} 


{-65.5} 


{-42.7} 


{5.776} 


{0.095} 






2.892 


77.3 


51.5 


-97.0 


67.7 


72.4 


-65.1 


-39.3 


5.382 


0.928 


1.324 


A2 


[3.779] 


[77.9] 


[51.3] 


[-98.8] 


[68.1] 


[72.8] 


[-65.5] 


[-38.3] 


[5.361] 


[1.189] 


[2.009] 




(3.715) 


(77.0) 


(51.6] 


(-97.4) 


(67.4) 


(72.8) 


(-65.4) 


(-38.6) 


(5.436) 


(0.655) 


(1.034) 




"3.730" 


"77.1" 


"51.6" 


"-96.6" 


"67.6" 


"72.4" 


"-64.9" 


"-38.4" 


"5.385" 


"1.000" 


"1.370" 




3.736 


77.4 


57.5 


-97.9 


67.4 


73.2 


-65.5 


-38.7 


5.438 


1.000 


1.018 




5.155 


118.5 


37.5 


-107.7 


71.0 


171.1 


-63.7 


57.6 


5.762 


1.202 






{5.004} 


{111.4} 


{37.6} 


{-101.9} 


{64.7} 


{164.9} 


{-61.1} 


{53.9} 


{5.586} 


{2.240} 






5.104 


118.6 


36.4 


-108.5 


69.9 


168.3 


-61.1 


53.4 


5.594 


2.303 


1.014 


A3 


[5.149] 


[121.2] 


[36.1] 


[-112.1] 


[71.9] 


[168.7] 


[-63.0] 


[55.9] 


[5.653] 


[0.585] 


[0.073] 




(5.079) 


(119.1) 


(36.8) 


(-110.2) 


(69.3) 


(168.7) 


(-60.9) 


(52.4) 


(5.702) 


(1.795) 


(0.553) 




"5.105" 


"120.0" 


"36.7" 


"-109.9" 


"71.0" 


"169.9" 


"-60.7" 


"52.5" 


"5.610" 


"2.864" 


"1.513" 




5.094 


119.1 


36.7 


-110.0 


69.6 


168.6 


-61.9 


54.8 


5.673 


2.864 


0.095 
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Table 2. Cont. 



Conformers 


R 


P a 


O a 


<Pi 


<P2 


q>3 


<P4 


<Ps 


t 1 


AE 


AG 




4.592 


116.2 


64.7 


-106.0 


73.2 


178.6 


65.0 


-56.8 


7.630 


1.726 






{4.487} 


{113.4} 


{64.9} 


{-100.7} 


{67.1} 


{178.0} 


{61.3} 


{-52.8} 


{8.026} 


{2.887} 






4.517 


117.5 


65.1 


-106.7 


71.7 


-179.5 


62.2 


-52.9 


7.390 


2.765 


1.457 


A4 


[4.552] 


[119.3 


[65.0] 


[-109.9] 


[73.7] 


[-179.3] 


[64.3] 


[-55.3] 


[7.309] 


[1.030] 


[0.432] 




(4.500) 


(118.5) 


(65.1) 


(-107.9) 


(71.0) 


(-179.6) 


(62.0) 


(-52.8) 


(7.480) 


(2.309) 


(1.009) 




"4.500" 


"117.5" 


"65.1" 


"-107" 


"71.0" 


"-179.6" 


"62.0" 


"-52.8" 


"7.480" 


"3.256" 


"1.871" 




4.515 


118.4 


65.1 


-108.2 


71.4 


-1 79. 7 


62.9 


-54.0 


7.467 


3.256 


0.503 


A5 


5.119 


123.3 


36.2 


-81.4 


-71.2 


-170.1 


60.7 


-53.6 


5.166 


2.952 


1.976 


A6 


4.839 


130.3 


50.2 


-93.9 


-178.3 


174.3 


-60.6 


54.4 


7.387 


3.870 


2.472 


A7 


4.342 


110.7 


69.0 


-100.5 


73.3 


-114.2 


-57.0 


53.3 


4.342 


4.088 


2.821 


A8 


5.314 


99.1 


61.8 


-108.8 


70.2 


174.4 


179.8 


123.1 


5.490 


4.529 


3.020 


A9 


4.316 


72.3 


29.9 


-100.9 


86.9 


94.8 


-69.7 


165.3 


8.401 


4.540 


3.284 


A10 


4.317 


98.8 


29.8 


-100.4 


87.1 


62.4 


118.4 


-119.6 


8.402 


4.541 


3.295 


All 


4.316 


98.8 


29.8 


-100.5 


87.1 


94.7 


-69.8 


165.6 


8.403 


4.541 


3.293 


A12 


5.248 


97.3 


61.6 


-107.1 


71.8 


176.2 


-179.4 


177.3 


6.700 


4.778 


3.132 


A13 


5.137 


115.6 


37.5 


-109.6 


69.9 


168.6 


-72.0 


169.7 


6.500 


4.912 


3.329 


A14 


3.804 


130.6 


83.0 


-137.1 


68.2 


-147.3 


63.6 


64.2 


7.916 


4.944 


4.017 


A15 


3.851 


140.4 


75.7 


-114.8 


-115.6 


140.7 


-61.6 


63.4 


7.669 


5.082 


3.917 


A16 


4.561 


120.1 


66.8 


-112.4 


70.5 


-178.3 


63.9 


60.5 


7.527 


5.103 


3.411 


A17 


4.400 


99.7 


34.4 


-112.0 


67.2 


57.8 


44.9 


48.3 


6.225 


5.139 


4.217 


A18 


5.204 


119.3 


39.3 


-110.7 


69.8 


173.8 


-64.3 


-63.3 


7.512 


5.216 


3.730 


A19 


3.411 


68.6 


64.3 


-92.1 


98.3 


-114.1 


64.4 


39.1 


7.760 


5.354 


4.890 


A20 


4.387 


136.6 


10.3 


-88.2 


158.8 


79.4 


-71.0 


172.6 


8.042 


5.378 


3.969 j 


A21 


3.562 


167.7 


60.8 


-94.2 


-126.1 


83.2 


-74.5 


178.0 


6.300 


5.384 


4.494 


A22 


4.486 


118.5 


67.8 


-112.8 


70.8 


-179.7 


71.7 


-171.9 


8.081 


5.432 


3.764 


A23 


5.209 


88.5 


51.3 


-104.5 


77.3 


91.3 


180.0 


-69.9 


6.218 


5.476 


3.817 


A24 


4.381 


115.3 


11.5 


-102.7 


-162.7 


-77.3 


73.0 


-170.7 


8.685 


5.514 


4.179 


A25 


4.948 


80.5 


59.5 


-102.5 


67.2 


74.6 


176.3 


-68.2 


6.425 


5.588 


4.277 


A26 


4.655 


80.0 


78.1 


-99.5 


74.6 


-104.0 


-175.2 


72.4 


7.000 


5.633 


4.093 
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Table 2. Cont. 



Conformers 


R 


P a 


O a 


<Pi 


<P2 


<P3 


<p 4 


<Ps 


t 1 


AE 


AG 


A27 


4.609 


83.1 


78.2 


-102.2 


71.5 


-107.6 


-174.8 


179.6 


5.473 


5.662 


4.371 


A28 


5.146 


128.0 


36.4 


-81.4 


-70.5 


-167.5 


71.8 


-170.7 


6.229 


5.675 


4.396 


A29 


4.653 


84.8 


51.8 


-103.7 


70.9 


-109.1 


-174.5 


-74.1 


6.731 


5.845 


4.325 


A30 


5.235 


90.1 


49.7 


-102.9 


80.0 


96.4 


179.4 


-173.8 


7.848 


5.850 


4.085 


A31 


4.400 


99.7 


34.4 


-112.0 


67.2 


57.8 


44.9 


48.3 


6.225 


5.890 


4.508 


A3 2 


5.220 


123.3 


38.2 


-81.4 


-70.2 


-172.6 


64.6 


64.8 


4.345 


5.909 


4.708 


A3 3 


4.415 


114.1 


12.6 


-104.1 


-163.1 


-74.6 


69.3 


69.4 


7.384 


5.917 


4.782 


A3 4 


5.955 


104.4 


31.3 


-95.4 


179.0 


-178.1 


179.2 


72.9 


6.080 


6.174 


4.503 


A3 5 


5.898 


124.0 


30.3 


-96.6 


179.3 


179.3 


179.9 


179.6 


6.676 


6.195 


4.544 


A3 6 


5.956 


124.3 


31.3 


-96.4 


179.5 


178.1 


-179.0 


-72.9 


7.716 


6.254 


4.762 


A3 7 


4.565 


102.4 


34.8 


-112.6 


66.4 


60.5 


56.3 


-178.8 


7.684 


6.290 


5.085 


A3 8 


4.273 


110.2 


72.9 


-105.8 


71.1 


-103.2 


-66.6 


167.8 


6.458 


6.432 


4.855 


A3 9 


4.970 


138.6 


46.2 


-100.0 


-179.8 


172.1 


-72.6 


166.4 


8.071 


6.610 


4.755 


A40 


5.011 


104.7 


46.2 


-97.2 


-179.7 


-176.1 


72.9 


-166.1 


7.023 


6.779 


5.097 


A41 


5.030 


104.4 


48.4 


-95.2 


179.6 


-177.9 


64.3 


61.5 


8.220 


7.100 


5.519 


A42 


5.722 


115.4 


19.7 


-99.3 


-176.9 


-84.1 


-179.2 


70.4 


6.405 


7.372 


5.289 


A43 


5.350 


87.6 


45.5 


-90.3 


138.5 


-88.9 


-177.6 


-70.3 


7.007 


7.483 


7.510 


A44 


5.386 


92.4 


40.2 


-91.6 


146.1 


-91.1 


178.5 


173.8 


7.133 


7.627 


5.567 


A45 


5.699 


118.8 


18.3 


-101.1 


-173.8 


-84.4 


-178.6 


176.1 


7.994 


7.665 


5.672 


A46 


5.765 


120.0 


17.5 


-101.1 


-170.2 


-80.2 


-173.9 


-69.9 


6.680 


7.684 


5.636 


A47 


4.832 


106.9 


33.0 


-107.0 


-164.9 


-72.0 


-62.7 


178.1 


8.554 


8.451 


6.636 




3.928 


73.4 


53.4 


75.1 


-146.7 


95.5 


-68.8 


72.8 


5.812 


0.072 






{3.885} 


{74.5} 


{50.3} 


{70.1} 


{-151.5} 


{94.2} 


{-68.1} 


{76.8} 


{6.204} 


{0.056} 






3.880 


73.4 


52.8 


73.5 


-145.4 


90.9 


-69.3 


76.7 


6.746 


0.083 


0.151 


Bl 


[3.927] 


[73.6] 


[53.1] 


[74.4] 


[-145.4] 


[92.6] 


[-69.8] 


[74.6] 


[5.485] 


[0.124] 


[0.822] 




(3.855) 


(73.5) 


(55.6) 


(74.2) 


(-145.7) 


(90.5) 


(-68.9) 


(76.7) 


(5.595) 


(0.061) 


(0.108) 




"3.870" 


"73.3" 


"51.1" 


"73.1" 


"-144.2" 


"89.6" 


"-70.0" 


"77.7" 


"5.374" 


"0.099" 


"0.195" 




3.882 


73.7 


52.7 


74.0 


-145.8 


91.5 


-69.3 


76.0 


5.601 


0.099 


0.090 
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Table 2. Cont. 



Conformers 


R 


P a 


O a 


<Pi 


<|>2 


<P3 


<|>4 


<P5 


t 1 


AE 


AG 




3.770 


77.2 


52.0 


97.0 


-66.7 


-74.1 


66.6 


36.1 


5.348 


1.216 






{3.726} 


{17.9} 


{50.4} 


{95.6} 


{-66.8} 


{-75.8} 


{65.5} 


{44.0} 


{5.664} 


{0}' 






2.905 


77.6 


51.2 


96.4 


-68.5 


-72.0 


65.1 


40.6 


5.340 


0.87 


1.251 


B2 


[3.784] 


[78.0] 


[50.8] 


[98.5] 


[-68.9] 


[-72.6] 


[65.4] 


[40.4] 


[5.321] 


[1.118] 


[1.900] 




(3.713) 


(77.3) 


(51.2) 


(96.8) 


(-68.2) 


(-72.4) 


(65.4) 


(39.9) 


(5.404) 


(0.612) 


(0.975) 




"3.727" 


"77.5" 


"51.2" 


"96.2" 


"-68.4" 


"-72.9" 


"64.0" 


"39.8" 


"5.250" 


"0.968" 


"1.316" 




3.736 


77.8 


51.1 


97.3 


-68.2 


-72.8 


65.4 


40.3 


5.398 


0.968 


0.941 




5.157 


117.3 


37.3 


106.1 


-71.7 


-171.1 


63.5 


-57.6 


5.951 


1.157 






{5.007} 


{110.9} 


{37.4} 


{100.4} 


{-65.2} 


{-165.1} 


{60.9} 


{-53.8} 


{5.998} 


{2.179} 






5.104 


117.9 


36.2 


107.1 


-70.3 


-168.1 


60.9 


-53.3 


5.839 


2.257 


0.903 


B3 


[5.149] 


[120.9] 


[35.9] 


[111.5] 


[-72.1] 


[-168.6] 


[62.9] 


[-55.9] 


[5.872] 


[0.542] 


[Of 




(5.080) 


(118.1) 


(36.5) 


(108.2) 


(-69.7) 


(-168.2) 


(60.8) 


(-53.3) 


(5.921) 


(1.748) 


(0.412) 




"5.106" 


"119.5" 


"36.5" 


"108.6" 


"-7 1.3" 


"-169.7" 


"60.6" 


"-52.4" 


"5.814" 


"2.819" 


"1.394" 




5.094 


118.6 


36.6 


108.9 


-69.9 


-168.4 


61.7 


-54.7 


5.907 


2.819 


0.001 




4.529 


95.5 


65.4 


78.3 


70.8 


-178.7 


63.6 


-54.4 


6.470 


2.031 






{3.673} 


{66.1} 


{78.5} 


{56.6} 


{63.8} 


{-147.2} 


{63.0} 


{-52.4} 


{6.649} 


{1.707} 






3.817 


68.7 


79.7 


56.7 


64.6 


-149.1 


63.8 


-50.2 


5.870 


2.722 


2.366 


B4 


[4.300] 


[85.9] 


[70.0] 


[69.9] 


[70.9] 


[-168.2] 


[63.2] 


[-52.1] 


[6.144] 


[1.331] 


[0.867] 




(3.814) 


(68.9) 


(79.4) 


(56.4) 


(64.5) 


(-150.0) 


(63.4) 


(-48.6) 


(5.913) 


(2.523) 


(2.080) 




"3.839" 


"68.8" 


"79.2" 


"55.9" 


"65.4" 


"-150.2" 


"63.6" 


"-48.6" 


"5.680" 


"3.179" 


"2.791" 




3.858 


70.3 


78.6 


58.0 


64.6 


-151.4 


63.4 


-49.6 


5.949 


3.179 


1.557 




4.598 


115.9 


64.6 


105.4 


-72.8 


-178.2 


-64.9 


56.7 


7.911 


1.693 






{4.507} 


{112.9} 


{64.3} 


{99.6} 


{-66.5} 


{-177.0} 


{-61.1} 


{52.8} 


{8.520} 


{2.832} 






4.526 


117.5 


64.8 


106.4 


-71.4 


180.0 


-62.0 


52.8 


7.735 


2.757 


1.463 


B5 


[4.559] 


[119.2] 


[64.9] 


[109.6] 


[-73.3] 


[179.7] 


[-64.2] 


[55.2] 


[7.629] 


[1.012] 


[0.420] 




(4.510) 


(118.4) 


(64.8) 


(107.6) 


(-70.7) 


-179.9 


(-61.8) 


(52.7) 


(7.826) 


(2.289) 


(0.978) 




"4.524" 


"117.6" 


"64.7" 


"106.6" 


"-72.3" 


"-179.9" 


"-6 1.7" 


"52.1" 


"7.517" 


"3.262" 


"1.920" 




4.524 


118.3 


64.8 


707.6 


-71.1 


-179.8 


-62.7 


53.8 


7.793 


3.262 


0.470 


B6 


5.119 


123.1 


36.3 


81.0 


70.9 


170.1 


-60.7 


54.1 


5.460 


2.993 


2.035 


B7 


4.991 


113.3 


29.3 


107.2 


-74.7 


-87.7 


-56.9 


46.5 


6.027 


3.636 


2.427 
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Table 2. Cont. 



Conformers 


R 


P a 


0 3 


<Pi 


<P2 


<P3 


<p 4 


<Ps 


t 1 


AE 


AG 


B8 


3.582 


165.9 


56.3 


82.7 


65.4 


73.0 


-75.2 


178.4 


6.463 


4.240 


3.708 


B9 


3.793 


70.4 


79.1 


57.2 


68.0 


-93.9 


-49.8 


60.8 


5.069 


4.263 


3.900 


BIO 


4.325 


99.2 


29.1 


98.3 


-88.0 


-94.7 


69.5 


-166.2 


8.268 


4.496 


3.334 


Bll 


5.318 


99.2 


61.7 


108.3 


-70.8 


-175.3 


179.6 


71.1 


5.640 


4.508 


3.006 


B12 


5.002 


132.0 


29.1 


77.5 


74.7 


88.4 


57.0 


-46.2 


5.222 


4.791 


3.824 


B13 


5.138 


115.0 


37.3 


108.2 


-70.2 


-168.4 


71.8 


-170.5 


6.133 


4.845 


3.212 


B14 


5.246 


98.4 


61.7 


108.6 


-71.5 


-176.6 


179.2 


177.1 


6.553 


4.852 


3.244 


B15 


5.312 


99.2 


62.1 


108.6 


-70.7 


-175.8 


-178.5 


-70.6 


7.964 


4.893 


3.413 


B16 


4.407 


99.6 


33.7 


110.2 


-68.1 


-57.9 


-44.6 


-48.7 


5.939 


5.042 


4.065 


B17 


5.205 


118.2 


39.0 


109.3 


-69.6 


-172.9 


64.6 


63.9 


7.429 


5.101 


3.551 


B18 


3.562 


167.3 


60.9 


93.8 


125.8 


-83.0 


74.6 


-178.0 


6.069 


5.456 


4.582 


B19 


5.277 


89.9 


49.5 


103.3 


-79.4 


-94.8 


-179.7 


69.2 


6.054 


5.577 


3.960 


B20 


4.876 


78.3 


61.1 


101.0 


-65.9 


-75.2 


179.3 


-69.0 


6.102 


5.694 


4.494 


B21 


5.253 


91.0 


49.4 


103.8 


-80.2 


-97.6 


-178.6 


174.8 


8.034 


5.950 


3.922 


B22 


4.534 


101.8 


35.0 


112.6 


-66.2 


-59.3 


-55.7 


179.6 


7.559 


6.169 


5.032 


B23 


3.146 


83.0 


95.1 


73.1 


67.3 


-116.9 


63.5 


64.3 


5.334 


6.305 


6.199 


B24 


4.610 


77.9 


77.5 


79.0 


72.8 


-102.2 


-173.9 


-175.2 


5.317 


6.617 


5.315 


B25 


4.365 


150.0 


33.5 


76.4 


67.0 


59.5 


40.5 


41.6 


5.882 


6.677 


5.427 


B26 


4.618 


108.0 


52.5 


77.2 


68.5 


-105.1 


-173.1 


-79.3 


5.308 


6.689 


5.500 


B27 


5.184 


143.6 


51.8 


78.9 


76.5 


89.5 


-178.9 


-69.4 


6.563 


6.718 


5.374 


B28 


4.635 


114.4 


78.3 


79.8 


73.8 


-101.1 


-174.8 


72.3 


7.031 


6.831 


5.483 


B29 


5.194 


141.1 


52.7 


77.5 


74.5 


87.3 


175.3 


68.0 


5.203 


6.997 


5.596 


B30 


5.129 


142.9 


52.3 


79.1 


76.5 


90.0 


-179.5 


-175.2 


7.031 


7.015 


5.609 


B31 


4.964 


136.4 


32.2 


77.2 


71.0 


78.8 


64.1 


-178.9 


5.033 


7.508 


6.051 



a Notation used from ref. [49]; b AE = 0 = -811.214731 a.u. at B3LYP/6-311++G(3df,pd) level; c AE = 0 = -810.723396 a.u. at B3LYP/6-31G(d,p) level; 
d AG = 0 = -810.764701 a.u.; e AE = 0 = -810.449499 a.u. at 03LYP/6-31G(d,p) level; f AE = 0 = -810.419555 a.u. at B3LYP/cc-pVDZ level; S AG = 0 = -810.460766 a.u.; 
h AE = 0 = -810.760402 a.u. at B3PW9 1/6-3 lG(d,p) level; ' AG = 0 = -810.801662 a.u.; j AE = 0 = -810.444324 a.u. at B972/6-31G(d,p) level; k AG = 0 = -810.485454 a.u.; 
1 AE = 0 = -808.643998 a.u. at MP2/6-31G(d) level; ra AG = 0 = -810.491798 a.u. 
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Table 3. Single point calculations at the MP2/6-31G(d,p)//B3LYP/6-31G(d,p) level. 



Conformer 



A£7kcalmol 



,-i 



Al 



0.143 



A2 



0.082 
2.553 
3.112 
0.122 



A3 



A4 



Bl 



B2 



0 



B3 
B4 



2.491 
1.962 



B5 



3.071 



The conformers differ in general very little in energy. Thus, in our calculations 78 optimized 
conformers were found within the electronic energy range AE = 0-8.5 kcal/mol, and Gibbs energy 
range AG = 0-6.6 kcal/mol with respect to the global minimum. This range of values of AG is smaller 
than that calculated in dG, in dT [50], 0-7 kcal/mol and in dU [31], 0-9 kcal/mol. 

Only four conformers are found within the electronic energy range AE = 0-1.0 kcal/mol (by 
criterium of AE + ZPE), Table 2, with cpi -727-97° as g~ and 70796° as g + by MP2. Among these 
conformers, Bl has the highest dipole moment 6.20 D although very close to that in Al, 6.18 D. These 
conformers are slightly favored in a polarizable environment with water. 

Another seven conformers appear within the electronic energy range AE = 1.0-3.0 kcal/mol, three 
values were anti (cpi ca. -100°) and four syn (cpi ca. 90°). The anti structures are the expected forms in the 
natural nucleosides that form the nucleotides and polynucleotides in biological systems [33,35,46^18]. 
The ratio anti/syn in conformers A is 0.8 in the low-energy group (<3 kcal/mol) but it increases up to 
1.7 in the 3.0-8.5 kcal/mol range. 

The global minimum by DFT methods corresponds to the conformer denoted as Al (Figure 4) and 
it appears stabilized by an intramolecular H-bond. The optimised bond lengths and natural NBO 
atomic charges on this conformer are collected in Figure 5. This global minimum by criterium of AE + 
ZPE agree well to that obtained by criterion of AG, but differs of that obtained by MP2/6-31G(d) and 
MP2/6-31G(d,p)//B3LYP/6-31G(d,p), Table 3 and Figure S3. It is because of the small difference in 
energy between both forms Al and B2. This global minimum in the syn form by MP2 is in accordance 
to that obtained in other nucleosides [33-35] but differs of the anti form expected for the natural 
nucleosides and nucleotides in biological systems [43]. The second most stable conformer is Bl with 
values of cpi = 70°, cp 2 = "152° and cp 3 = 94° by MP2. 

Figure 4 shows the six best optimum conformers selected in the two ranges of (pi: three are A 
(A1-A3), and three are B (B1-B3). The values of the intramolecular H-bonds and the most important 
structural angles of each conformer are also included. Figures 6-9 show the distribution of the 78 
optimised conformers according to their energies, exocyclic torsional angles cpi-cps, the angles P, O, 
and the vector R. The best significant conformers are pointed in these figures. 
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Figure 4. Geometry of the six most optimum conformers selected for each rotation angle 
(pi determined in Al conformer of Nl tautomer of acyclovir at B3LYP/6-31G(d,p) level. 
The values of the strongest intramolecular H-bonds are also included. 




Figure 5. Natural atomic charges and optimum bond lengths in conformer B2 of tautomer 
Nl of acyclovir at B3LYP/6-31G(d,p) and MP2/6-31G(d,p) levels, in fuchsia and green 
colours, respectively. 



Geometry NBO charges 
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Figure 6. Distribution of the 78 optimum stable calculated conformers at the B3LYP/6- 
31G(d,p) level in tautomer Nl of acyclovir according to their exocyclic torsional angle cpi 
and their: (a) relative electronic energy AE + ZPE correction; (b) relative Gibbs energy 
AG; and (c) the distance R. The most stable conformers of each type are pointed. 
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Figure 7. Distribution of the 78 optimum stable calculated conformers in tautomer Nl of 
acyclovir, according to the values of the exocyclic torsional angles: 92. 93- 94 and 95 versus 
the angle 91. 
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Figure 8. Distribution of the 78 optimum stable calculated conformers in tautomer Nl of 
acyclovir according to the values of the angles <J> and [3, and the dipole moment \i versus the 
angle (pi. 
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Figure 9. 3D plots with the relative energies of the 78 optimum stable conformers 
according to the values of the exocyclic torsional angles: cpi. cp 2 - 93 and cp 4 versus relative 
electronic energy AE + ZPE. 




3.3. Conformational Angle Analysis 

An overall examination of the five exocyclic torsional angles and two bond angles, defining the 
conformational space in tautomer Nl of ACV, leads to conclude the following, Figures 7-9: 

(i) The dihedral angle cpi presents a bimodal distribution: (-137° < -74°) 47 conformers A (anti) 
and (57° < 113°) 31 conformers B (syn). Anti forms prevail in number over syn ones and cover 
a slight wider range of cpi values. This anti orientation has been found in the crystal of several 
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nucleosides [51-53], and it is the form for biological activity [54]. Although, in many 
nucleoside analogues syn and anti forms have similar energy, however, the global minimum 
corresponds in general to the syn form, as in AZT [35]. In purine nucleosides there is relatively 
little restraint to rotation about the N9-C1' bond [55]. In ACV the high population of the anti 
forms and the wider range of (pi values are factors that facilitate its antiviral activity. The sugar 
ring of the natural nucleoside dG can be mimics by the similar exocyclic angles of ACV. 

(ii) The torsional cp 2 angle has a clear bimodal distribution: 65° < (p 2 < 180° (49 conformers, 34 are 
A and 15 are B) and -180° < q> 2 < "66° (29 conformers, 13 are A and 16 are B), Figure 7. 

(iii) The (p 3 angle also has a bimodal conformation: 58° < (p 3 < 180° (32 conformers, 22 are A and 
10 are B) and -180° < q> 3 < "58° (46 conformers, 25 are A and 21 are B). 

(iv) The cp4 angle has a trimodal distribution: -173° < (p 4 < 175° (28 conformers, 16 are A and 12 
are B), 41° < q> 4 < 118° (28 conformers, 17 are A and 11 are B), and -75° < q> 4 < -45° 
(22 conformers, 14 are A and 8 are B). 

(v) The torsional cp 5 angle also has a trimodal distribution: -166° < cps < 165° (29 conformers, 
19 A and 10 B), 39° < cp 5 < 123° range (29 conformers, 16 A and 13 B), and -120° < cp 5 < "46° 
(20 conformers, 12 A and 8 B). 

(vi) In the bond angle [3 has been observed a regular distribution in the 69° < P < 168° range, Figure 
8. This large range indicates the high flexibility of the chain with multiple 05'H orientations. 
The most stable conformers with AE < 1 kcal-moi -1 have the [3 angle with values between 72° 
and 78°, while the conformers with higher AE have some values larger than 140°. 

(vii) The bond O angle also shows a regular distribution in the 10° < O < 83° range for conformers 
A, and 29° < O < 95° for conformers B. The value of this angle appears between 51° and 54° 
for the four best conformers with AE < 1 kcal/moi -1 . In general, the most stable conformers do 
not have large values of O, neither low values. 

Finally, the R vector has a large range of values, between 2.905A and 5.956 A, and it determines 
the distance of the OH end of the chain from the base plane. 

3.4. Guanine Moiety 

In the analysis of the six most stable conformers, the base heterocycle has a very small 
non-planarity, in general with torsional angles lower than 1°, with the exception of those involving N3 
and N9 atoms with a value higher than 2°. It is due to the strong intramolecular H5'--N3 hydrogen 
bond which produces a ring deformation around N3. Thus, among the dihedral angles describing the 
base heterocycle nonplanarity, we have defined: v 0 (N1-C2=N3-C4), Vi (C2=N3-C4=C5), v 2 
(N3-C4=C5-C6), v 3 (N7=C8-N9-C4) and v 4 (C6-N1-C2=N3). 

An analysis of the six most stable conformers (Figure 4) shows that the highest base heterocycle 
nonplanarity corresponds to conformer A2 with the torsional angles Vi (2.8°), v 2 (-2.3°) and Vo 
(-1.3°). The other conformers show deviations higher than 1° in the following torsional angles: 
Conformer Al in v 2 (-1.2°) and v 3 (-1.1°), conformer A3 in Vi (1.5°) and v 2 (-1.3°), conformer Bl in 
v 0 (1.3°) and Vi (1.1°), and conformer B2 in v 4 (1.2°). 
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Comparing the bond lengths of ACV (conformer Al) with those of dG, it is observed that ACV has 
several bonds of the guanine moiety larger than the corresponding bonds in dG, while in the remaining 
bonds are smaller. The largest difference appears in C4-N3. 

3.5. The Side Chain 

The most important characteristic of the structure of ACV is the conformation of the side chain that 
is attached to N9, which is characterized by the exocyclic torsional angles cpi to cps, by the bond angles 
P and O that determine the position of the chain from the base, and by the distance R. The bonds along 
this chain can be in trans respect to each other and thus it gives an extended zig-zag structure or, by 
contrast, one or several alternate bonds can be in gauche in such a way that it resembles at least 
partially a portion of the furanose moiety of dG [49,56]. In this last case, it is capable of adopting 
conformations resembling a portion of the pentose ring, a factor which undoubtedly plays an important 
role in their biological activities [57]. Thus, three types of structures appear: two that have the side 
chain partially folded (conformers A and B), and one structure which has all bonds in trans orientation 
and leading to an almost planar zig-zag arrangement [49]. 

3. 6. Intramolecular H-Bonds 

Several authors have studied the intramolecular H-bonds in related nucleosides, in special using 
AIM method [58,59]. However, the absence of the furanose ring in ACV reduces largely the number 
of possible H-bonds. Therefore, only two intramolecular H-bonds may be observed in the main 
conformers of ACV: (i) The hydroxyl hydrogen H5'(05') and the guanine's position 3 nitrogen atom, 
H5'---N3, and (ii) The hydroxyl hydrogen H5'(05') and the guanine's position 2 oxygen atom, 
H5'-02'. 

Figure 4 shows these types of H-bonds in the six most stable conformers. H-bond (i) appears in 
conformers Al, Bl, B2 and A2; while (ii) is observed in conformers B3 and A3 with values of 
2.346 A and 2.350 A, respectively. H-bond (i) is stronger than (ii), and it gives a great stability to the 
structure. Thus, the conformers with the H-bond (i) are the most stable, AE < 1 kcal-mol , while the 
conformers with the H-bond (ii) have higher energy, AE > 2 kcal-mol 1 . 

3.7. Natural NBO Atomic Charges 

The calculated values in conformer Al appear collected in Figure 5. The largest negative charge 
corresponds to N2 and 05' atoms, -0.90 and -0.84 e, respectively by MP2 (where e is the charge of an 
electron). The next atom with large negative charge is Nl (ca. -0.7 e). The value of this charge is 
slightly higher than in dG, -0.69 e. The main effect of the NH2 group is a remarkable increment in the 
positive charge on C2, 0.81 e, because of the high negative charge on the amino nitrogen, -0.90 e. 
Consequently, a noticeable increase of the negative charge on N3 is observed, -0.73 e (-0.70 e in dG). 
The electron-rich sites of the guanine moiety are involved in H-bonds, with N3 and 02' acting as 
single acceptors. 

In the nitrogen N9 the negative charge (-0.52 e) is lower than on Nl, N3 and N7, but slightly 
higher than in dG. It is because in ACV the bonding to the chain increases the negative charge on N9. 
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The value of the charge on N7 is important because in the neutral form of anti-tumour platinum drugs, 
the platinum atom has a strong preference for nitrogen N7 rather than for oxygen atoms of the base for 
its coordination [51]. Also, the N7 position in DNA is the most open to attack. Nl position is also 
important because when deprotonation of the weakly acid ACV occurs, the metal binding site changes 
to Nl, which is the formally deprotonated site [60]. Raman spectra of related nucleosides in H2O 
indicate that the site of deprotonation in basic solutions is Nl, while the site of proto nation in acidic 
solutions is N3 and N7, the same sites that in its phosphorylated form. These results are useful for 
identification and characterization of its structure in natural occurring biopolymers [61]. 

The negative charge on the oxygen atoms is high, -0.65 e on 02' and -0.69 e on 06. The values in 
02' and 05' are slightly higher than in dG, by B3LYP -0.603 e and -0.761 e, respectively. By contrast, 
in 06 the negative charge is lower than in dG, -0.592 e. 

C6 is the atom with the highest positive charge, 0.829 e by MP2, in concordance to the high 
negative charge on 06, i.e., it is the most reactive. With a slightly lower value appears C2 (0.812 e), 
and with much lower values H5'(05), H1(N1), H2(N2) andH2'(N2). The remaining hydrogen atoms have 
much less positive charge and they are less reactive. 

3.5. Solid State Simulation 

The structure of ACV (CSD code CEHTAK) in the solid state was determined by X-ray 
diffraction [53,60-63]. Birnbaum et al. [49,62] obtained the crystals that belong to monoclinic space 
with three independent molecules (A, B and C), together with two water molecules [49]. Molecules A 
and B have similar conformation while that of molecule C is different, Table 4. In molecules A and B 
the angle (pi and 92 are in the preferred g~ form, but in molecule C 92 is in trans. It means that molecules A 
and B showed a partially folded conformation of the side chain, while molecule C appears extended [62]. 
Molecules A and B correspond to conformer A5 in the isolated state, while molecule C is represented by 
conformer A34. Molecules A and B are ca. 3 kcal-mol more stable than molecule C. The difference in 
energy between molecule A and B is small and it is due to the flexibility of the side chains. 

Table 4. A comparison of the most important structural parameters in the dimer and 
tetramer forms calculated at the B3LYP/6-31G(d,p) level in acyclovir molecule with those 
in the crystal. The torsional angles are in degrees and R in A. 
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We have simulated this arrangement through a dimer and tetramer forms, Figure 10. The calculated 
value of the inter- and intramolecular H-bonds distances and the total energy is included in the Figure. 
The calculated bond lengths, bond angles and torsional angles are collected in Tables SI and S2. 
The simulated structure is tightly bonded but not so tight as it the reported in the crystal. There are four 
donor sites: H1(N1), H2(N2), H2'(N2) and H5', and five possible acceptors: 06, N3, N7, 02' and 05'. 



Figure 10. The optimized dimer and tetramer forms in conformer Al of acyclovir at the 
B3LYP/6-31G(d,p) level. The H-bonds observed are in A and the E (RB+HF-LYP) in a. u. 




Dimer -1621.895988 a.u. Tetramer -3243.852071 a.u. 



The dimer, and tetramer forms were built with the most stable conformer Al. In the dimer there are 
two intermolecular H-bonds between the two molecules: Nl {in molecule A} and 06 {in molecule B} 1 .733 A, 
and 06{in A} and Nl {in B} 1.732 A. The tetramer shows seven H-bonds: (i) two between molecule A 

and B: H(Ni){in molecule A} and 06{in molecule B} 1.709 A, and 06{in molecule A} and H(Ni){in molecule B} 
1.847 A; (ii) There is only a H-bond between molecule B and C: H (0 5 ){in B} and N3{in C}, 1.945 A; 
(iii) There are two H-bonds between molecule C and D: H(Ni){in molecule C} and 06{in molecule D}, 
1.741 A, and 06{inC} and H( N i){inD}, 1.792 A. (iv) At list there are two H-bonds between molecule D 
andB: N7{in molecule D} and H^^in B}, 1.972 A. 

Comparing the calculated bond lengths and angles of molecule A in the dimer with the X-ray 
data [62], Tables S2 and S3, it is observed the good agreement in the values of our simulated structure. 
The greater differences appear in C4=C5 (0.025 A), C2-N2 (0.023 A), C3'-C4' and N1-C6 (0.016 A), 
and the difference is almost null in C6=0, C5-N7, N7=C8 and C4-N9. A slight larger error is obtained 
when the comparison theory-experiment is carried out with the molecule B of the dimer. This error is 
slight reduced when the comparison is carried out with the simulated molecule A of the tetramer. In 
the angles the difference is almost null, with similar values of all the angles. Only several angles in the 
side chain show some noticeable differences, i.e., C3'-C4'-05', 4.7°. 

Other four anhydrous forms of ACV and a new hydrate have been characterized by X-ray 
diffraction, with significant differences in the intermolecular H-bonding networks among the ACV 
forms [63-65]. A study of the solution forms of ACV shows that ACV can exist as polymorphic and 
pseudopolymorphic solvates [66]. We have simulated other dimer, trimer and pentamer forms with 
different conformers and collected in Figure S4. 
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4. Conclusions 

In the present work we have shown a comprehensive compendium of the possible conformers in 
tautomer Nl of ACV. The geometries and values of the properties presented here appear to the most 
accurate to date. The most important findings of the present work are the following: 

(1) Five tautomers of ACV were identified and fully optimized. At room temperature only 
tautomer Nl (48.1%), OHC (37.7%) and OHT (14.3%) have a noticeable population. It is very 
small in tautomers N3 and N7, less than 0.05%. 

(2) The relative energies of the five tautomers appear related to the dipole moment and to the 
torsional angle cpi. The least stable tautomer N3 in the isolated state has the highest jj, and thus, 
it is the most favoured in a polarisable environment with water. 

(3) In the isolated state the most stable tautomer is Nl by both B3LYP and MP2 methods. In this 
tautomer, and through a rotation of cpi, q>2, q>3, (p4 and cp5 exocyclic torsional angles, 78 
optimized stable conformers were identified, two syn and two anti failing into the 0-1 kcal-moF 1 
AE + ZPE energy range. 

(4) The calculated most stable conformer by all DFT levels corresponds to Al, while by MP2 is 
B2. In the nitrogen atoms and in the 02' and 05' oxygen atoms of conformer Al appear a 
higher reactivity than in the corresponding natural nucleoside deoxyguanosine. 

(5) The distribution of all the conformers according to the ranges of stability of the characteristic 
torsional angles was established. The values obtained indicate the flexible nature of ACV, 
which is higher than dG. An increase of the stability appears when the side chain is near to the 
purine base, with a value of R that fails into 3.925-2.892 A range, and an angle O close to 54°. 

(6) Only two intramolecular H-bonds may be observed in the main conformers of ACV, in contrast 
to the six H-bond types calculated in dG. It leads to a flexibility higher in ACV than in dG. 

(7) The solid state was simulated through a dimer and tetramer forms. An excellent agreement with 
the X-ray crystal data was obtained, which indicates the good accuracy of the theoretical 
methods used. 

Author Contributions 

This research corresponds to the PhD Thesis of M.C. Alvarez-Ros, in which M. Alcolea Palafox is 
the director. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 



1. Schloemer, G.C.; Han, Y.-K.; Harrington, P.J. Preparation of acyclovir using 1,3-dioxolane. 
U.S. Patent 5567816, 22 October 1996. 



Pharmaceuticals 2014, 7 



718 



2. Lazrek, H.B.; Baddi, L.; Smietana, M.; Vasseur, J. -J.; Sebti, S.; Zahouily, M. One-pot synthesis of 
antiviral acyclovir and other nucleosides derivatives using doped natural phosphate as Lewis acid 
catalyst. Nucleos. Nucleot. Nucleic Acids 2008, 27, 1 107-1 1 12. 

3. Gao, H.; Mitra, A.K. Synthesis of acyclovir, ganciclovir and their prodrugs: A review. Synthesis 
2000, 3, 329-351. 

4. De Clercq, E. Antiviral drugs: Current state of the art. J. Clinical Virol. 2001, 22, 73-89. 

5. Emmert, D.H.; Lancaster, M.D. Treatment of common cutaneous herpes simplex virus infections. 
Am. Fam. Physician. 2000, 61, 1697-1607. 

6. Fuchs, J.; Celum, C; Wang, J.; Hughes, J.; Sanchez, J.; Cowan, F.; Reid, S.; Delany-Moretlwe, S.; 
Corey, L.; Wald, A. Clinical and virologic efficacy of herpes simplex virus type 2 suppression by 
acyclovir in a multicontinent clinical trial. J. Infect Dis. 2010, 201, 1 164-1 168. 

7. Piret, J.; Boivin, G. Resistance of herpes simplex viruses to nucleoside analogues: Mechanisms, 
prevalence, and management. Antimicrob. Agents Chemother. 2011, 55, 459-472. 

8. Harnden, M.R.; Wyatt, P.G.; Boyd, M.R.; Sutton, D. Synthesis and antiviral activity of 
9-alcoxypurines. 1.9-(3-hydroxypropoxy)- and 9-[3-hydroxy-2-(hydroxymethyl)propoxy]purines. 
J. Med. Chem. 1990, 33, 187-196. 

9. Golankiewicz, B.; Ostrowski, T.; Goslinski, T.; Januszczyk, P.; Zeidler, J.; Baranowski, D.; 
de Clercq, E. Fluorescent tricyclic analogues of acyclovir and ganciclovir. A structure-antiviral 
activity study. J. Med. Chem. 2001, 44, 4284-4287. 

10. Lisco, A.; Vanpouille, C; Tchesnokov, E.P.; Grivel, J.C.; Biancotto, A.; Brichacek, B.; Elliott, J.; 
Fromentin, E.; Shattock, R.; Anton, P.; et al. Acyclovir is activated into a HIV-1 reverse 
transcriptase inhibitor in herpes virus-infected human tissues. Cell Host Microbe. 2008, 4, 260-270. 

11. Vanpouille, C; Lisco, A.; Derudas, M.; Saba, E.; Grive, J.-C; Brichacek, B.; Scrimieri, F.; 
Schinazi, R.; Schols, D.; McGuigan, C; et al. A new class of dual-targeted antivirals: 
Monophosphorylated acyclovir prodrug derivatives suppress both human immunodeficiency virus 
type 1 and herpes simplex virus type 2. J. Infect Dis. 2010, 201, 635-643. 

12. Tanton, C; Weiss, H.A.; LeGoff, J.; Changalucha, J.; Clayton, T.C.; Ross, DA.; Belec, L.; 
Hayes, R.J.; Watson- Jones, D. Patterns of herpes simplex virus shedding over 1 month and the 
impact of acyclovir and HIV in HSV-2-seropositive women in Tanzania. Sex Transmit. Infect. 
2011,57, 406-411. 

13. Delany, S.; Mlaba, N.; Clayton, T.; Akpomiemie, G; Capovilla, A.; Legoff, J.; Belec, L.; Stevens, W.; 
Rees, H.; Mayaud, P. Impact of acyclovir on genital and plasma HIV-1 RNA in HSV-2/HIV-1 
co-infected women: a randomized placebo-controlled trial in South Africa. AIDS. 2009, 23, 461^169. 

14. Prichard, M.N.; Kern, E.R.; Hartline, C.B.; Lanier, E.R.; Quenelle, D.C. CMX001 Potentiates the 
efficacy of acyclovir in herpes simplex virus infections. Antimicrob. Agents Chemother. 2011, 55, 
4728-4734. 

15. Kalachandra, S.; Takamata, T.; Lin, D.M.; Snyder, E.A.; Webster-Cyriaque, J. Stability and release 
of antiviral drugs from ethylene vinyl acetate (EVA) copolymer. J. Mater. Sci: Mater. Med. 2006, 
17, 1227-1236. 

16. Gavira, J.M.; de la Fuente, M.; Navarro, R.; Hernanz, A. Normal coordinate analysis of 
acycloguanosine. J. Molec. Struct. 1997, 410-411, 425-429. 



Pharmaceuticals 2014, 7 



719 



17. Birnbaum, G.I.; Johansson, N.G.; Shugar, D. Conformations of acyclonucleosides: Crystal 
structure of 9-(4-hydroxybutyl)guanine, an analogue of acyclovir. Nucleos. Nucleot. Nucleic Acids 
1987, 6, 775-783. 

18. Brandi-Blanco, M.P.; Choquesillo-Lazarte, D.; Dominguez-Martin, A.; Gonzalez-Perez, J.M.; 
Castineiras, A.; Niclos-Gutierrez, J. Metal ion binding patterns of acyclovir: Molecular recognition 
between this antiviral agent and copper(II) chelates with iminodiacetate or glycylglycinate. 
J. Inorg. Biochem. 2011, 105, 616-623. 

19. Barcelo-Oliver, M.; Terron, A.; Garcia-Raso, A.; Fiol, J.J.; Molins, E.; Miravitlles, C. Ternary 
complexes metal [Co(II), Ni(II), Cu(II) and Zn(II)]-ortho-iodohippurate (I-hip)-acyclovir. X-ray 
characterization of isostructural [(Co, Ni or Zn)(I-hip) 2 (ACV)(H 2 0)3] with stacking as a recognition 
factor. J. Inorg. Biochem. 2004, 98, 1703-1711. 

20. Brovarets, O.O. ; Yurenko, Y.P.; Dubey, I.Y.; Hovorun, D.M. Can DNA-binding proteins of 
replisome tautomerize nucleotide bases? Ab initio model study. J. Biomol. Struct. Dyn. 2012, 29, 
1101-1109. 

21. Brovarets, O.O.; Hovorun, D.M. Prototropic tautomerism and basic molecular principles of 
hypoxanthine mutagenicity: An exhaustive quantum-chemical analysis. J. Biomol. Struct. Dyn. 
2013, 31, 913-936. 

22. Plass, M.; Kristl, A.; Abraham, M.H. Spectroscopic investigation of tautomeric equilibria in the 
guanine derivatives of acyclovir. J. Chem. Soc. Perkin Trans II 1999, 264, 26-41. 

23. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. 
Phys. Rev. A 1988, 38, 3098-3100. 

24. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 
1993, 98, 5648-5652. 

25. Handy, N.C.; Cohen, A.J. Left-right correlation energy. Mol. Phys. 2001, 99, 403-412. 

26. Hoe, W.-M.; Cohen, A.J.; Handy, N.C. Assessment of a new local exchange functional OPTX. 
Chem. Phys. Lett. 2001, 341, 319-328. 

27. Xu, X.; Zhang, Q.; Muller, R.P.; Goddard III, WA. An extended hybrid density functional 
(X3LYP) with improved descriptions of nonbond interactions and thermodynamic properties of 
molecular systems. J. Chem. Phys. 2005, 122, 14105-14114. 

28. Lee, C; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into 
a functional of the electron density. Phys. Rev. 1988, B37, 785-789. 

29. Perdew, J.P.; Wang, Y. Accurate and simple analytic representation of electron-gas correlation 
energy. Phys. Rev. 1992,545, 13244-13249. 

30. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; 
Montgomery, J.A., Jr.; Vreven, T.; Kudin, K.N.; Burant, J.C.; et al. Gaussian 03, Revision B.04. 
Gaussian Inc.: Pittsburgh, PA, USA, 2003. 

31. Yurenko, Y.P.; Zhurakivsky, R.O.; Ghomi, M.; Samijlenko, S.P.; Hovorun, D.M. Ab Initio 
comprehensive conformational analysis of 2'-deoxyuridine, the biologically significant DNA minor 
nucleoside, and reconstruction of its low- temperature matrix Infrared spectrum. J. Phys. Chem. B 
2008, 112, 1240-1250. 



Pharmaceuticals 2014, 7 



720 



32. Yurenko, Y.P.; Zhurakivsky, R.O.; Ghomi, M.; Samijlenko, S.P.; Hovorun, D.M. Comprehensive 
conformational analysis of the nucleoside analogue 2'-[3-deoxy-6-azacytidine by DFT and MP2 
calculations. J. Phys. Chem. B 2007, 111, 6263-6271. 

33. Alcolea Palafox, M.; Iza, N.; de la Fuente, M.; Navarro, R. Simulation of the first hydration shell 
of nucleosides D4T and Thymidine: structures obtained using MP2 and DFT methods. J. Phys. 
Chem. B 2009, 113, 2458-2476. 

34. Alcolea Palafox, M.; Iza, N. Structure of the antiviral stavudine using quantum chemical methods: 
complete conformational space analysis, 3D potential energy surfaces and solid state simulation. 
J. Molec. Struct. 2012, 1028, 181-195. 

35. Alcolea Palafox, M.; Talaya, J. Hydration analysis of antiviral agent AZT by means of DFT and 
MP2 calculations. J. Phys. Chem. B 2010, 114, 15199-15211. 

36. Alcolea Palafox, M.; Posada-Moreno, P.; Villarino-Marin, A.L.; Martinez-Rincon, C; 
Ortuno-Soriano, I.; Zaragoza-Garcia, I. DFT Calculation of four new potential agents muscarinic 
of bispyridinium type: structure, synthesis, biological activity, hydration, and relations with the 
potents W84 and DUO-30. J. Comput. Aided Molec. Des. 2011, 25, 145-161. 

37. Alcolea Palafox, M. Anticancer drug IUdR and other 5-halogen derivatives of 2'-deoxyuridine: 
conformers, hydrates and structure-activity relationships. Struct. Chem. 2014, 25, 53-69. 

38. Alcolea Palafox, M.; Nielsen, O.F.; Lang, K.; Garg, P.; Rastogi, V.K. Geometry and vibrational 
spectra of 5-substituted uracils. Asian Chem. Letts. 2004, 8, 81-93. 

39. Alcolea Palafox, M.; Rastogi, V.K. Quantum Chemical predictions of the vibrational spectra of 
polyatomic molecules. The uracil molecule and two derivatives. Spectrochim. Acta. 2002, 58, 
411-440. 

40. Alcolea Palafox, M. The Prediction of Vibrational Spectra: The Use of Scale Factors. In Recent 
Research Developments in Physical Chemistry; Transworld Research Network: Trivandrum, 
India, 1998; Volume 2. 

41. Alcolea Palafox, M. Scaling factors for the prediction of vibrational spectra. I. Benzene molecule. 
Int. J. Quantum Chem. 2000, 77, 661-684. 

42. Alcolea Palafox, M.; Iza, N.; Gil, M. The hydration effect on the uracil frequencies: an 
experimental and quantum chemical study. J. Molec. Struct. 2002, 585, 69-92. 

43. Saenger, W. Principles in Nucleic Acid Structure; Springer Verlag: New York, NY, USA, 1984. 

44. Carpenter, J.E.; Weinhold, F. Analysis of the geometry of the hydroxymethyl radical by the "different 
hybrids for different spins" natural bond orbital procedure. J. Molec. Struct. 1988, 169, 41-62. 

45. Wesolowski, T.A.; Godzik, A.; Geller, M. Calculations of the conformational properties of 
Acyclonucleosides. Part I. Stable conformations of acyclovir. Acta Biochim. Polonica. 1987, 34, 
111-122. 

46. Tamara Molina, A.; Alcolea Palafox, M. Structure and conformational analysis of the anti-HIV 
AZT 5'-aminocarbonylphosphonate prodrug using DFT methods. Chem. Phys. 2011, 387, 1 1-24. 

47. Alcolea Palafox, M. Molecular structure differences between the antiviral Nucleoside Analogue 
5-iodo-2'-deoxyuridine and the natural nucleoside 2'-deoxythymidine using MP2 and DFT 
methods: Conformational analysis, crystal simulations, DNA pairs and possible behaviour. J. 
Biomol. Struct. Dyn. 2014, 32, 831-851. 



Pharmaceuticals 2014, 7 



721 



48. Alcolea Palafox, M.; Iza, N. Tautomerism of the Natural Thymidine Nucleoside and in the 
Antiviral Analogue D4T. Structure and influence of an Aqueous Environment Using MP2 and 
DFT Methods. Phys. Chem. Chem. Phys. 2010, 12, 881-893. 

49. Birnbaum, G.I.; Cygler, M.; Kusmierek, J.T.; Shugar, D., Structure and conformation of the potent 
antiherpes 9-(2-hydroxyethoxymethyl)guanine (acycloguanosine). Biochem. Biophys. Res. Comm. 
1981, 103, 968-974. 

50. Yurenko, Y.P.; Zhurakivsky, R.O.; Ghomi, M.; Samijlenko, S.P.; Hovorun, D.M. How many 
conformers determine the thymidine low-temperature matrix infrared spectrum? DFT and MP2 
quantum chemical study. J. Phys. Chem. B. 2007, 111, 9655-9663. 

51. Parthasarathy, R.; Kim, H. Conformation and sandwiching of bases by azido groups in the crystal 
structure of 3'-azido-3'-deoxy-thymidine (AZT), an antiviral agent that inhibits HIV reverse 
transcriptase. Biochem. Biophys. Res. Comm. 1988, 152, 351-358. 

52. Camerman, A.; Mastropaolo, D.; Camerman, N. Azido thymidine: crystal structure and possible 
functional role of the azido group. Proc. Natl. Acad. Sci. USA 1987, 84, 8239-8242. 

53. Dyer, I.; Low, J.N.; Tollin, P.; Wilson, H.R.; Howie, RA. Structure of 3'-deoxythymidine, AZT. 
Acta Crystallog. 1988, C44, 767-769. 

54. Painter, G.R.; Andrews, C.W.; Furman, PA. Conformation and local environment of nucleotides 
bound to HIV type 1 reverse transcriptase (HIV-1 RT) in the ground state. Nucleos. Nucleot. 
Nucleic Acids. 2000, 19, 13-29. 

55. Thewalt, U.; Bugg, C.E.; Marsh, R.E. The crystal structure of guanosine dihydrate and inosine 
dihydrate. Acta Cryst. 1970, B26, 1089-1101. 

56. Kazimierczuk, Z.; Stolarski, R.; Shugar, D. Acyclonucleosides: acyclobenzimidazole nucleoside 
and nucleotide analogues and conformations of the acyclic chains by means of NMR spectroscopy. 
Acta Biochim. Polonica. 1984, 31, 33-48. 

57. Birnbaum, G.I.; Stolarski, R.; Kazimierczuk, Z.; Shugar, D. Solid state and solution conformations 
of l-(P-D-2'-3'-secoribofuranosyl)-5,6-dichlorobenzimidazole, an acyclonucleoside analogue. 
Can. J. Chem. 1985, 63, 1215-1221. 

58. Fidanza, N.G.; Sosa, G.L.; Lobayan, R.M.; Peruchena, N.M. Topological analysis of the 
electronic charge density in nucleoside analogues derivatives of the AZT. Effects of X-H---0 and 
X-H---F intramolecular H-bonds. J. Molec. Struct. 2005, 722, 65-78. 

59. Yurenko, Y.P.; Zhurakivsky, R.O.; Samijlenko, S.P.; Hovorun, D.M. Intramolecular CH...0 
hydrogen bonds in the Al and Bl DNA-like conformers of canonical nucleosides and their 
Watson-Crick pairs. Quantum chemical and AIM analysis. J. Biomol. Struct. Dyn. 2011, 29, 51-65. 

60. Kliifers, P.; Mayer, P. The sodium salt of Nl-deprotonated inosine, Na .C10H11N4O5 .2.5H 2 0. 
Acta Crystallogr. 1996, C52, 2970-2972. 

61. Medeiros, G.C.; Thomas, G.-J., Jr. Raman studies of nucleic acids IV. Vibrational spectra and 
associative interactions of aqueos inosine derivatives. Biochim. Biophys. Acta 1971, 247, 449-462. 

62. Birnbaum, G.I.; Cygler, M.; Shugar, D. Conformational features of acyclonucleosides: structure 
of acyclovir, an antiherpes agent. Can. J. Chem. 1984, 62, 2646-2652. 

63. Lutker, K.M.; Quinones, R.; Xu, J.; Ramamoorthy, A.; Matzger, A.J. Polymorphs and hydrates of 
acyclovir. J. Pharm. Sci. 2011, 100, 949-963. 



Pharmaceuticals 2014, 7 



722 



64. Sohn, Y.T.; Kim, S.H. Polymorphism and pseudopolymorphism of acyclovir. Arch. Pharm. Res. 
2008, 37,231-234. 

65. Barboza, F.; Vecchia, D.D.; Tagliari, M.P.; Silva, M.A.S.; Stulzer, H.K. Differential scanning 
calorimetry as a screening technique in compatibility studies of acyclovir extended release 
formulations. Pharm. Chem. J. 2009, 43, 363-368. 

66. Kristl, A.; Srcic, S.; Vrecer, F.; Sustar, B.; Vojnovic, D. Polymorphism and pseudopolymorphism: 
Influencing the dissolution properties of the guanine derivative acyclovir. Int. J. Pharm. 1996, 
139, 231-235. 



© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



